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We investigate the energy band structure of the superfluid flow of ultracold dilute Fermi gases in 
a one-dimensional optical lattice along the BCS to BEC crossover within a mean-field approach. In 
each side of the crossover region, a loop structure (swallowtail) appears in the Bloch energy band 
of the superfluid above a critical value of the interaction strength. The width of the swallowtail is 
largest near unitarity. Across the critical value of the interaction strength, the profiles of density 
and pairing field change more drastically in the BCS side than in the BEC side. It is found that 
along with the appearance of the swallowtail, there exists a narrow band in the quasiparticle energy 
spectrum close to the chemical potential and the incompressibility of the Fermi gas consequently 
experiences a profound dip in the BCS side, unlike in the BEC side. 
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Ultracold atoms in optical lattices attract much inter- 
est because the controllability of both the lattice geome- 
try and the interatomic interaction is such that they serve 
as testing beds for various models [1]. For Bose-Einstein 
condensates (BECs), it has been pointed out that the 
interaction can change the Bloch band structure drasti- 
cally, causing the appearance of a loop structure called 
"swallowtail" in the energy dispersion [2, 3]. This is due 
to the competition between the external periodic poten- 
tial and the nonlinear mean-field interaction: the former 
favors a sinusoidal band structure, while the latter tends 
to make the density smoother and the energy dispersion 
quadratic. When the nonlinearity wins, the effect of the 
external potential is screened and a swallowtail energy 
loop appears [4]. This nonlinear effect requires the exis- 
tence of an order parameter and, consequently, the emer- 
gence of swallowtails can be viewed as a peculiar mani- 
festation of superfluidity in periodic potentials. 

The problem of swallowtails can be even more impor- 
tant in Fermi superfluids due to the possible wide impli- 
cations for various systems in condensed matter physics 
and nuclear physics. Indeed, extensive recent research 
activities are devoted to the simulation of solid states us- 
ing cold Fermi gases and the behavior of cold fermions in 
optical lattices can lead to interesting analogies with su- 
perconductors and superconductor superlattices. In ad- 
dition, our work may have implications also for superfluid 
neutrons in neutron stars, especially those in "pasta" 
phases (see, e.g., Ref. [5] and references therein) in neu- 
tron star crusts [6] , where nuclei form a crystalline lattice 
in which superfluid neutrons can flow. However, unlike 
the Bose case [2-4, 7-9], little has been studied in this 
problem so far and a fundamental question of whether 
or not swallowtails exist along the crossover from the 
Bardeen-Cooper-SchriefFer (BCS) to BEC states is still 
open. In this context, our work is aimed at showing the 



existence and the conditions for emergence of swallow- 
tails in Fermi superfluids and presenting the unique fea- 
tures which make them different from those in bosons. 

We consider a two-component unpolarized dilute Fermi 
gas made of atoms of mass m interacting with s-wave 
scattering length and subject to a one-dimensional 
(ID) optical lattice of the form 14;xt(r) = sE'r sin^ q-QZ = 
Vo sin^ qbz. Here, Vq = sEs, is the lattice height, s is the 
lattice intensity in dimensionlcss units, i^R = h?q^/2m 
is the recoil energy, qb = n/d is the Bragg wave vector 
and d is the lattice constant. We compute the energy 
band structure of the system by solving the Bogoliubov- 
de Gennes (BdG) equations [10]: 
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where H{r) — -h?V^/2m. + \4xt(i") - Ui{r) and 
Vi{T) are quasiparticle wavcfunctions, and is the cor- 
responding quasiparticle energy. The chemical poten- 
tial ^ is determined from the constraint on the parti- 
cle number N ~ 2^./ \vi{r)\'^ dr and the pairing field 
A(r) should satisfy a self-consistency condition A(r) = 
—gJ2i where g is the coupling constant for the 

s-wave contact interaction which needs to be renormal- 
ized. In the presence of a supercurrent with wavevec- 
tor Q ~ P/h {\P\ < Podge = ^9b/2) moving in the 
z-direction [11], one can write the quasiparticle wave- 
functions in the Bloch form as Ui{r) = Ui{z)e'^^^e'^^'^ 
and Vi{v) = Vi{z)e~'^^^e^^''^ leading to the pairing field 
as A(r) — e*^'5^A(z). Here A{z), Ui{z), and Vi{z) are 
complex functions with period d and the wave vector 
{\kz\ < ^b) lies in the first Brillouin zone. This Bloch de- 
composition transforms Eq. (1) into the following BdG 
equations for Ui(z) and Ui{z) : 
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FIG. 1: (Color online) (a) Energy E per particle as a function 
of the quasimomentum P for various values of 1/fcpas; (b) 
half-width of the swallowtails along the BCS-BEC crossover. 
These results are obtained for s = 0.1 and Ef/Eb, = 2.5. 
The quasimomentum Pedge = %b/2 fixes the edge of the first 
Brillouin zone. The dotted line in (b) is the half-width in 
a EEC obtained by solving the GP equation; it vanishes at 
l/kpas ~ 10.6. 



where 



HQiz) 



2m 



Here, k\ = k"^ + ky and the label i represents the wave 
vector k as well as the band index. 

In the following, we mainly present the result for s = 
Vq/E-r = 0.1 and E-p/En = 2.5 as an example, where 
E-p = fi^/cp/(2m) and kp = (STr^np)^/^ are the Fermi 
energy and momentum of a uniform free Fermi gas of 
density ng. These values fall in the range of parameters 
of feasible experiments [12]. 

We first compute the energy per particle in the low- 
est Bloch band as a function of the quasimomentum P 
for various values of l/kpag- The results in Fig. 1(a) 
show that the swallowtails appear above a critical value 
of l/fepfts where the interaction energy is strong enough 
to dominate the lattice potential. In Fig. 1(b), the half- 
width of the swallowtails from the BCS to the BEC side is 
shown. It reaches a maximum near unitarity (l/fepOs = 
0). In the far BCS and BEC limits, the width vanishes 
because the system is very weakly interacting and the 
band structure tends to be sinusoidal. When approach- 
ing unitarity from either side, the interaction energy in- 
creases and can dominate over the periodic potential, 
which means that the system behaves more like a trans- 
lationally invariant superfluid and the band structure fol- 
lows a quadratic dispersion terminating at a maximum 
P larger than Pcdgo- In the BEC side, we compare the 
results of our BdG calculations with those of the Gross- 
Pitaveskii (GP) equation for bosons of mass mi, = 2m 
interacting with scattering length ab = 2as [13] in an op- 
tical lattice 2Kxt(z): -(?iV4m)V2$(r) 2T4xt(z)*(r) + 
(87r/i^as/2TO)|<i>(r)p$(r) = ^B$(r), where $(r) is a sin- 
gle macroscopic wavefunction describing the BEC. As it 
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FIG. 2: (Color online) The parameter region {Ep/E-n, 
l/kpas) where the swallowtails appear for s = 0.1. Sym- 
bols on the vertical line at Ep/E^^ — 2.5 correspond to the 
cases shown in Fig. 1(a). 



should be, the difference between the two curves becomes 
vanishingly small in the deep BEC regime. With our set 
of parameters (s = 0.1 and E-p/Eji = 2.5), the width of 
the swallowtail predicted by the GP equation vanishes at 
l/kpas ~ 10.6. 

Whether the swallowtail appears or not depends on 
three parameters: the interaction parameter l/fcpas, the 
lattice intensity s, and the ratio between the Fermi 
energy and the recoil energy Ep/En- In Fig. 2, we 
fix s = 0.1 and show the parameter region {Ep/E-n, 
l/Zcpfls) where we find swallowtails in the BCS side of 
the crossover {l/kpUg < 0) [14]. We see that for weaker 
interaction (i.e., larger values of I/ZcfIosI) higher den- 
sities (larger values of Ep) are required to create swal- 
lowtails, as expected. The s-dependence of critical val- 
ues of l/kpas in the BCS side is much weaker than the 
1/s scaling behavior in the BEC side. In the BCS side, 
as far as we have checked in 0.1 < s < 0.5, the criti- 
cal value of l/kpGs changes within 30% at Ep/E^ =- 5 
and the change gets smaller with increasing Ep/Ep{. 
This weak dependence of the swallowtail region on s in 
the BCS side is due to the Fermi statistics: Provided 
Ep/Vo ~ {Ep/Epi)/ s is sufficiently larger than unity, the 
flow of the BCS condensate formed from fcrmions near 
the Fermi surface is not very sensitive to the presence of 
the lattice potential, while the flow of the BEC formed 
from bosonic dimers all at the bottom of the energy levels 
is more sensitive to it. 

Both the pairing field and the density exhibit interest- 
ing features in the range of parameters where the swal- 
lowtails appear. This is particularly evident at the Bril- 
louin zone boundary, P = Pedge- In Fig. 3, we show the 
magnitude of the pairing field |A(z)| and the density n[z) 
calculated at the minimum [z = 0) and at the maximum 
(z = ±d/2) of the lattice potential. In general, n(z) and 
|A(z)| take maximum (minimum) values where the exter- 
nal potential takes its minimum (maximum) values (for 
the full profiles, see Supplemental Material [15]). The fig- 
ure shows that |A(d/2)| remains zero in the BCS regime 
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FIG. 3: (Color online) Profiles of (a) the pairing field |A(z) 
and (b) the density n{z) along the change of l/fcpas for P = 
Podge in the case of s = 0.1 and Ef/E-e{, = 2.5. The values 
of |A(2)| and n{z) at the minimum {z — 0, blue □) and at 
the maximum {z = ±d/2, red x) of the lattice potential are 
shown. The vertical dotted lines show the critical value of 
l/kpUs above which the swallowtail exists. The dotted curve 
in (a) shows |A[ in the uniform system. 



until the swallowtail appears at l/fepag w —0.62. Then it 
increases abruptly to values comparable to |A(0)|, which 
means that the pairing field becomes almost uniform at 
P = Pedge in the presence of swallowtails. As regards 
the density, we find that the amplitude of the density 
variation, n(0) — n(d/2), exhibits a pronounced maxi- 
mum near the critical value of l/fepOs- In contrast, in 
the BEC side, the order parameter and the density are 
smooth monotonic functions of the interaction strength 
even in the region where the swallowtail appears. At 
P = Podge, the solution of the GP equation for bosonic 
dimers gives the densities ni,{0) ~ moil + Vo/2ni,oUo) 
and nb{d/2) = nbo{l - Vo/2ribo[/o) with Vo/2nboUo = 
(37r/4)(s£'R/£'F)(l/fcFas), where nto is the average den- 
sity of bosons and C/q = 'inh'^ab/'mt [2, 16, 17]. Near 
the critical value of l/k-pQs, unlike the BCS side, the 
nonuniformity just decreases all the way even after the 
swallowtail appears. The local density at z ~ d/2 is zero 
until the swallowtail appears in the BEC side while it 
is nonzero in the BCS side irrespective of the existence 
of the swallowtail. The qualitative behavior of |A(z)| 
around the critical point of l/fepfls is similar to that of 
nt,{z) because nb{r) = (m^as/87r)|A(r)|^ [18] in the BEC 
limit. 

The quasiparticle energy spectrum plays an important 
role in determining the properties of the Fermi gas. Here 
we show that the emergence of swallowtails in the BCS 
side and for Ep / E'r ^ 1 is associated with peculiar struc- 
tures of the quasiparticle energy spectrum around the 
chemical potential. In the presence of a superflow moving 
in the z direction with wavevector Q, the quasiparticle 
energies are given by the eigenvalues in Eq. (2). Since 
the potential is shallow (s ^ 1), some qualitative results 
can be obtained even ignoring Vcxtiz) except for its pe- 
riodicity. With this assumption we obtain 
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FIG. 4: (Color online) (a) Lowest three Bloch bands of the 
quasiparticle energy spectrum at fcx = for P = Pedge and 
l/fcras = —0.62. Thin black dashed lines labeled by Ts 
show the approximate energy bands obtained from Eq. (3) 
using ^ ~ 2.66Sr and |A| ~ |A(0)| ~ 0.54£;r. (b) Incom- 
pressibility at P = Pedge around the critical value of 
l/fcras ~ —0.62 where the swallowtail starts to appear. The 
quantity is the incompressibility of the homogeneous free 
Fermi gas of the same average density. In both panels, we 
have used the values s = 0.1 and Ep/Eb. — 2.5. 



with I being integers for the band index. If Q = 
0, the I = band has the energy spectrum 
^y[{k^ + k1)/2m — nY + |A]2 which has a local maxi- 
mum at kz = k± ~ 0. When Q — Podgo/^, the spectrum 
is tilted and the local maximum moves to k^ ~ qB/2 pro- 
vided JA] ^ Ep (and Ep/E^ > 1). In the absence of the 
swallowtail, the full BdG calculation indeed gives a local 
maximum at kz = qB/2 [see Fig. 4(a), where the green 
dashed line in the region —1 < k^/qB < —0.5 and the red 
solid line in the region —0.5 < k^/qB < 1 correspond to 
the I = band] and the quasiparticle spectrum is sym- 
metric about this point, which reflects that the current 
is zero. As Ep/En increases, the band becomes flatter 
as a function of kz and narrower in energy. 

In Fig. 4(a), we show the quasiparticle energy spec- 
trum at k± = for P = Podge- When the swallowtail is 
on the edge of appearing, the top of the narrow band just 
touches the chemical potential /i [see the dotted ellipse 
in Fig. 4(a)]. Suppose l/kpUg is slightly larger than the 
critical value, so that the top of the band is slightly above 
fi. In this situation a small change of the quasimomen- 
tum P causes a change of /i. In fact, when P is increased 
from P = Podgo to larger values, the band is tilted and 
the top of the band moves upwards; the chemical po- 
tential n should also increase to compensate for the loss 
of states available. This implies dfi/dP > 0. (See also 
Supplemental Material [15].) On the other hand, since 
the system is periodic, the existence of a branch of sta- 
tionary states with dfi/dP > at P = Podgo implies the 
existence of another symmetric branch with dfi/dP < 
at the same point, thus suggesting the occurrence of a 
swallowtail structure. 

A direct consequence of the existence of a narrow 
band in the quasiparticle spectrum near the chemical 
potential is a strong reduction of the incompressibility 
= ndfi(n)/dn close to the critical value of l/Zcpfls in 
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the region where swallowtails exist in the BCS side [see 
Fig. 4(b)]. The dip of occurs in the situation where 
the top of the narrow band is just above /Lt for P = Pcdgc 
{l/kpag is slightly above the critical value). An increase 
of the density n has little effect on /x in this case, because 
the density of states is large in this range of energy and 
the new particles can easily adjust themselves near the 
top of the band by a small increase of fx. This implies 
that dfi{n)/dn is small and the incompressibility has a 
pronoimccd dip [19]. It is worth noting that in the BEC 
side the appearance of the swallowtail is not associated 
with any significant change of incompressibility. In fact, 
the exact solution of the GP equation gives = ni,QUo 
near the critical conditions for the occurrence of swal- 
lowtails, being a smooth and monotonic function of the 
interaction strength. 

Swallowtails may produce observable effects in the be- 
havior of Bloch oscillations [3, 20]. Since Bloch oscilla- 
tions have various important applications, such as preci- 
sion measurements of forces [21] and controlling the mo- 
tion of a wave packet [22, 23], a better understanding of 
swallowtails in bosonic and fermionic gases is certainly 
useful in these contexts. One may also exploit the tun- 
ability of the interaction and the peculiar dynamics of the 
superfluid in the lattice for different applications. For in- 
stance, by periodically sweeping a magnetic field across 
the critical region for the appearance of swallowtails, one 
can produce a time modulation of the shape of the lowest 
energy band between a sinusoidal form and a quadratic- 
like form. Since the absolute value of the group velocity 
dp[E{P) /N] of the latter is always larger than that of the 
former [see Fig. 1(a)], one could realize a directed motion 
of the gas by synchronizing the period of this modulation 
with the period of the Bloch oscillations (see also Supple- 
mental Material [15]). This may be experimentally more 
accessible in the BCS side than in the BEC side because 
the critical value of I/^fIosI is of order 1 for a wide range 
of s and E-p/E^. This new method would complement 
other proposals for realizing directed motion of atomic 
wave packets in ID optical lattices [23, 24]. 

In summary, we have predicted the existence of swal- 
lowtails in the energy band of superfluid fermions in a 
lattice and have pointed out some key features which 
make these swallowtails different from those in a BEC. 
The results arc obtained within a range of parameters 
compatible with current experiments [12]. We hope our 
predictions stimulate experiments aimed to observe swal- 
lowtails with Fermi gases. 
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Full profiles of the pairing field |A(2:) 
density n{z) 



and the interaction strengths are short, just enough, and suffi- 

ciently large for developing the swallowtails, respectively. 



In Fig. 3 of the paper, the values of |A(z)| and n{z) at 
the minimum and at the maximum of the lattice potential 
are shown. The full profiles of |A(z)| and n{z) along 
the lattice vector (z-direction) are given in the following 
Fig. 5. By increasing the interaction parameter l/fcpas, 
we find that the order parameter |A| at the maximum 
{z = ±d/2) of the lattice potential exhibits a transition 
from zero to nonzero values at the critical value of 1/k-pas 
at which the swallowtail appears [see also Fig. 3(a) of the 
paper]. Note that here we plot the absolute value of A; 
the order parameter A behaves smoothly and changes 
sign across zero. 
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FIG. 5: Profiles of the pairing field |A(2:)| and the density 
n{z) at l/kpaa = —0.8 (blue dotted), —0.6 (green dashed), 
and —0.4 (red solid) for P = Pcdgc in the case of s = 0.1 and 
Ef/Eb. = 2.5. The swallowtail starts to appear at a critical 
value of l/fcpas ~ —0.62. 

The blue dotted, green dashed, and red solid 
lines correspond to the cases where the interparticle 



Behavior of the Bloch Band near the value of the 
chemical potential with the change of the 
quasimomentum 

In Fig. 4(a) of the paper, the lowest three Bloch bands 
of the quasiparticle energy spectrum for P = Podge at 
k± = and l/fepOs = —0.62 are given. To visualize the 
behavior of the second band (the red solid line in that 
figure) near the value of the chemical potential with the 
change of the quasimomentum P, we show the cases of 
P/Pedge = (black dotted), 0.5 (green dashed), and 1 
(red solid) in the following Fig. 6. Notice that the sharp 
minima in the curves for P/Pedgo = and —0.5 arc due 
to avoided crossings with other bands. 
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FIG. 6: Bloch band of the quasiparticle energy spectrum 
around the chemical potential for P/Pedge = (black dotted), 
0.5 (green dashed), and 1 (red solid), respectively, at fc^ = 
and l/feptis ~ —0.62 in the case of s = 0.1 and Ef/Er = 2.5. 
The each horizontal line denotes the value of the chemical 
potential for the corresponding value of P. 



Application of the swallowtail band structure to a 
Directed motion 



In Fig. 7, we show a schematic diagram for our proposal 
of the directed motion of an atomic wave packet. The ba- 
sic idea is making use of the difference in the group veloc- 
ity VG between the quadratic-like band with a swallowtail 
and the sinusoidal-like band. By modulating l/kpUg in 
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FIG. 7: Schematic diagram of the group velocity vq during 
one period (Tmod) of modulation of the magnetic field (ac- 
cordingly, l/fcpaa). The green (blue) dashed line shows the 
group velocity along the time when l/feptts is set to the value 
for a quadratic (sinusoidal) band structure as depicted in Fig. 
1(a) of the paper. When l/kpas is switched between the two 
values in a proper way, the group velocity of a wave packet 
will follow the black solid line. In real situation, the black 
line around t = Tmod/S is a smooth curve rather than a sharp 
drop due to a continuous variation of l/fcpas connecting the 
two values. 



such a way that the dispersion is quadratic when vq is 
positive and sinusoidal when vg is negative, the net dis- 
placement over one period of modulation of l/k-pUs is 
positive and a directed motion can be produced. 



